Quantum optomechanics of a Bose-Einstein Antiferromagnet 
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We investigate the cavity optomechanical properties of an antiferromagnetic Bose-Einstein con- 
densate, where the role of the mechanical element is played by spin-wave excitations. We show 
how this system can be described by a single rotor that can be prepared deep in the quantum 
regime under realizable experimental conditions. This system provides a bottom-up realization of 
dispersive rotational optomechanics, and opens the door to the direct observation of quantum spin 
fluctuations. 
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PACS numbers: 42.50.Pq, 03.75.Mn 

Cavity optomechanics (COM) has witnessed consider- 
able advances in recent years, promising significant con- 
tributions in several topics of broad interest including 
the production and manipulation of macroscopic quan- 
tum states, the investigation of the crossover between 
quantum and classical regimes in systems of increasing 
size, and the realization of a wide range of applications 
from high-precision metrology to quantum information 
science and motion-mediated nonlinear photonics [D-(jl|. 

The first such systems featured optomechanical cou- 
plings with a linear dependence on the displacement of 
a vibrating mechanical element such as a mirror 
More recently, systems featuring quadratic optomechan- 
ical coupling have also been realized 0]. This advance 
was quite important since, in principle, quadratic op- 
tomechanical coupling allows the non-destructive mea- 
surement of the phonon number state of the oscillator, as 
well as observation of the discrete quantum jumps that 
accompany changes in phonon number 0,0]. In contrast, 
these measurements are impossible with only linear op- 
tomechanical coupling. 

A parallel development in COM was the realization of 
COM systems formed by ultracold atomic gases trapped 
inside driven high-Q resonators [H[ , the role of the macro- 
scopic mechanical element now being played by a col- 
lective excitation of the atoms. These systems can be 
thought of as "bottom-up" realizations of COM, since in 
practice one can readily prepare a degenerate atomic gas 
in its quantum ground state, and thus the motion of the 
fictitious mechanical element is initially in the quantum 
regime [HJ. This is in contrast to the more conventional 
"top-down" COM systems described above where despite 
considerable progress it remains in general an ongoing 
challenge to cool the optomechanically-coupled mechan- 
ical mode into the quantum regime. 

This Letter considers a bottom-up situation where the 
role of the mechanical element is played by spin-wave ex- 
citations of a macroscopic mode of an antiferromagnetic 
Bose-Einstein condensate (BEC), whose dynamics can be 
mapped to a quantum rotor model 0-0] ■ Although there 
have been multiple investigations of the cavity dynamics 



of a spin ensemble I3~14|, as well as an atom-chip- based 
realization of quadratic-coupling COM this is to our 
knowledge the first study of bottom-up COM where the 
atomic gas has spin degrees of freedom and the cavity 
confinement leads to the ability to resolve the quantum 
regime of an antiferromagnetic spin gas. 



Especially pertinent to the present study is Ref. [12j . 
which provides the first link that we know of between 
the fields of spin systems and COM control in the liter- 
ature. That paper studies the dynamics of a large quan- 
tum spin in an optical resonator and finds spin bistability 
and squeezing of light by establishing an analogy with 
a harmonic oscillator coupled linearly to a single cavity 
mode. In contrast, the present Letter describes an exper- 
imentally realizable process for the dramatic cooling of a 
collective spin-wave excitation mode by establishing an 
analogy with a torsional oscillator coupled quadratically 
to a single cavity mode. As already mentioned the form 
of the optomechanical coupling is crucial in determining 
which phenomena can be observed in such a system. 

We consider an atomic spin-1 Bosc condensate char- 
acterized by antiferromagnetic spin-exchange collisions, 
and described by the Hamiltonian (h = 1) 



(1) 



Here C2 is the antiferromagnetic spin coupling, N is 
the total particle number, and F — tpjFijipj is the to- 
tal spin operator, where are the spin-1 matrices 
and tpi (i = ±, 0) are the bosonic annihilation opera- 
tors for each spin component. The additional coupling, 
Q = (/ i B-B) 2 /(4Ahf), is the quadratic Zeeman shift, with 
external magnetic field B and hypcrfinc splitting A^f. 

Recently, it has been shown that there is an exact map- 
ping of the Hamiltonian (|T]) to a Hamiltonian describing 
a single quantum rotor [8|, |l5| , 



q 2 N 



V(0) = q{N + 3/2)8^0 + ^— sin 2 (20), (2) 
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where the conjugate operators 9 and L are the ro- 
tor's effective angular displacement and angular momen- 
tum, respectively, and the effective moment-of-inertia is 
/ = N/c 2 . An effective position space basis is speci- 
fied by defining the eigenstates of the Cartesian compo- 
nents \X) = (|+) - \-))/y/2, \Y) = (|+) + |-»/»V2, 
and \Z) = |0), where {|±), |0)} is the single-particle spin 
eigenbasis. 

We focus on the weak magnetic field limit q <C c 2 where 
the dynamics of the quantum rotor are dominated by the 
first term of V(9), and are thus localized around the poles 
of the Bloch sphere. By expanding V{9) to quadratic 
order about the minimum 9 = we find that the quantum 
rotor model of Eq. ([2]) is further reduced to one that 
describes a one-dimensional harmonic torsional oscillator 
characterized by 9 and L z , where [9, L z ] = i 0, H|, and 
9 is the polar angle of the total spin state. In that limit 
the approximate potential is 
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for which the ground state wave-function is 
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One finds however that the harmonic torsional oscilla- 
tor description is only valid when 9 <C 1. Combining 
this requirement with the initial weak-field condition, the 
parameter region where harmonic approximation of the 
quantum rotor is valid is 



1 < c 2 /q < 2A 2 



(5) 



When this condition is met, any finite field will localize 
the state of the spinor condensate about 9 = regard- 
less of its initial state [H [13] • In this so-called "Joseph- 
son regime" @ the rotor model can be used to predict 
several purely quantum mechanical phenomena such as 
small-number effects (c 2 /q = 20, A = 200 or 9 ~ 0.02) 
0, observation of spin collapse and revival (9 = 0.1) 
and ncmatic-mixing dynamics [l6j]. These effects are 
difficult to observe experimentally because for coupling 
strengths characteristic of antiferromagnetic BEC, mean- 
field effects will dominate unless one uses relatively small 
condensates, A ~ tens to hundreds of atoms Q. This 
requirement is in conflict with the requirement of rel- 
atively large condensates needed to characterize them 
through the usual probe of absorption imaging. As a 
result most studies of spinor BECs have focused on mean- 
field physics 

However, when an antiferromagnetic BEC is trapped 
in a driven high-Q optical cavity, the atom-photon cou- 
pling is greatly enhanced, and can thus be used as a sen- 
sitive experimental probe for the BEC. It follows that, 
under presently available experimental conditions (see, 



for example, Ref. [H|), one can measure the quantum 
dynamics of small condensates, and as we will show, even 
observe the quantum regime in large condensates. More 
specifically, such cavity confinement allows one to enter 
the quantum regime where an antiferromagnetic spin-1 
BEC is described by the quantum rotor model. 

When trapped inside a unidirectional ring cavity the 
antiferromagnetic BEC is described by the Hamiltonian 



H — H-a + He + Hl 
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where = L 2 /{2I) + {\/2)Iu 2 9 2 is the quantum rotor 
Hamiltonian, Hq describes the intracavity field and its 
interaction with the rotor, H^ = — ini,{a — a') describes 
the laser pumping with amplitude k l , and Ho e and H-y 
account for mechanical and input noise, respectively. 

The effective optomechanical coupling arises within 
Hq- In the rotating frame of the pump laser 

H c =(uj c - u L )a ] a + U (ip\.ip + + ^_ip-)a?a , (7) 

where wl and wc are the frequencies of the pump laser 
and the empty cavity resonance, respectively. The second 
term of Eq. ([7|) describes the far off-resonant coupling of 
the cavity field to the BEC, where U$ = g 2 / (cjl — w a ) is 
the single-photon light shift with dipole coupling g and 
atomic transition frequency ui & . Assuming that the cavity 
field is 7r-polarized, the atoms in the F = 1 ground-state 
manifold couple to the F' = 1 excited-state manifold. 
The selection rules for the corresponding transitions are 
Amp = 0, with the exception that the \F = l,mp = 
0) — > \F' = l,mp> = 0) transition is forbidden, which 
means that atoms in the spin-0 state do not couple to 
the cavity field . 

With the normalization condition ip+ t l-'+ + ip—i/*— + 
V'oV'o = A, the mapping —qipQipo — > V(9), and in the 
harmonic approximation Hq can be reexpressed as 



H c ~{-A + t B 9 2 )a)a, 



(8) 



where A = — (wc + U$N) is the static part of the 
pump-cavity detuning. The dynamic part of the pump- 
cavity detuning is readily identified as the quadratic op- 
tomechanical coupling term Hqm — £,e9 2 a 1f a, with 
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The most salient feature of Hom is the quadratic de- 
pendence on 9, which, in principle, allows one to perform 
quantum non-demolition measurements of the rotational 
energy of the oscillator as well as to observe the dis- 
crete quantum jumps that accompany a change in roton 
number [3, H| , thereby providing a unique probe of the 
quantum state of an antiferromagmetic BEC. We remark 
that, as in other studies of systems featuring quadratic 
optomechanical coupling 0, [j, one can adjust the 
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system parameters, in this case tune the magnetic field, 
to realize higher order couplings. For example 



H gi = H 6 2 - P9 4 , /3~(q- U (Ja)N/3. 



(10) 



where H S 2 is the quadratic coupling Hamiltonian dis- 
cussed above. 

We now turn to the system dynamics of Hg2 , which 
we study by considering the corresponding Heisenberg- 
Langcvin equations 
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Here 7 is the cavity damping rate; a m is the operator 
describing the input noise, characterized by a zero mean 
and Markovian correlations (a m (i)a m, t(i')) = S(t — t'), 
and e g n accounts for noise from a thermal reservoir and/or 
incoherent atomic scattering, with (s e n ) = and 
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where Dg is the intrinsic damping constant for the quan- 
tum rotor and fce is Boltzmann's constant. 

In the mean-field approximation, the steady-state so- 
lutions of the Heisenberg-Langevin equations are 
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An important feature of these solutions is the lack of 
multistability, a common hallmark of many cavity op- 
todynamical systems, including both the linear optoro- 
tational coupling case fl9| and the general cavity-spin 
coupling case 10h13|. The lack of a multistablc region 
plays a pivotal role in the enhanced cooling and trapping 
of the effective quantum rotor: it allows one to achieve 
stronger cooling and trapping by increasing the pump 
laser power without having to worry about instabilities 
arising due to the emergence of additional stable states. 
Of course, when the harmonic motion condition of Eq. ([5]) 
is severely violated, such as in the limit of high magnetic 
fields (q > C2) treated in Ref. 13|, strong bistability tends 
to occur. 

Similarly to the case of a linearly coupled rotating mir- 
ror (l9j . the quadratic optomechanical coupling of the 
quantum rotor results in an increased effective trapping 
frequency. We show this by replacing each system op- 
erator with its first-order fluctuation expansion 0(t) = 
O s + SO(t), writing the cavity field fluctuation in terms 
of its quadrature operators SXi = (5a + 5a^)/y/2, SX 2 = 
(5a — 5a^)/i^/2, and then linearizing the Heisenberg- 
Langevin equations with respect to fluctuations to obtain 

i>{t)=Rv(t) + A(t), v(t) = {S9,SL Z ,SX 1 ,6X 2 ) T , 
A(t) = (0, ef, V^/X[ n , y^X™), (13) 
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and X[ n = (a in + a in ^)/y/2 and Xf = (a in - a in >t)/i\/2 
are the quadrature operators for the input-field fluctua- 
tions. It is important to note that a steady state solution 
is only stable if the real part of each of the correspond- 
ing eigenvalues of R is non-positive, which can be easily 
confirmed by using the Routh-Hurwitz criterion (20j . 

After Fourier transform and subsequent algebraic ma- 
nipulation one solves for the first-order fluctuations 

(15) 

, v r , V27{(7-^l n 2H±AA»M} 
5X h2 [uj] = f — '■ . (16 
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Equation (|15[) is in the form of a linear response 

89[u] = xM4 n M> 

where the noise operator e™ plays the role of the external 
perturbation, and with 
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where x(w) is the susceptibility of a fictitious mechanical 
rotor characterized by the damping factor r c ff = Dg and 
the resonant frequency cj c ff = r/ujg, where we have defined 
the enhancement factor 



(18) 



In contrast to the linear-coupling case [la ], rj is indepen- 
dent of frequency, maximized at A = 0, and can be as 
large as 10 3 or even higher for typical parameter values. 
Stiffcr trapping and even more efficient cooling can be ex- 
pected for small values of q, a key factor for the detection 
of quantum fluctuation effects [7H9( . 

The effect of this enhancement on the mean roton num- 
ber can be found from the first-order correlation func- 
tions of the system. From Eq. (fTTj) we find 

d , o2 > 1 
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FIG. 1: Average roton occupation number n as a function 
of scaled static detuning for T= 2 /iK and 500 pK, Uo = 
2tt x 100 Hz, 7 = 2tt x 50 kHz, n L = 2tt x 3 MHz. 



regime of the rotor model, which has been previously 
unattainable. Furthermore, since the system we describe 
is within the reach of current experimental capabilities, 
our results can be of immediate use to researchers study- 
ing such diverse problems as spin-based nano-mechanical 
devices [23[ , squeezed rotors [l8| , hybrid spin- mirror en- 
tanglement, and high-precision control of spinor atoms, 
polar molecules, or even rotating BECs [24J. 
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Note added: After this manuscript was submitted, an- 
other paper 23j presented a related theoretical analysis 
on a torsional model of a cavity-spin system. 



where the term featuring n = [exp(u!g j ksT) — 1] , the 
mean thermal excitation number at temperature T, ac- 
counts for a shift of the first-order correlations due to con- 
tact with a thermal reservoir. The steady state solution 
to Eq. (fH)]) is stable and exhibits non-vanishing expecta- 
tions (9z)ss and {L 2 Z ) SS - The steady-state energy of the 
quantized rotor is written Eq — (L 2 Z ) SS /2I + Iojg(6^) ss /2. 
Using the steady-state solution of Eq. (TlT)|) one finds 
Eq = (n+ l/2)u' e , where u' g = (w 8 /2)[l + r?" 2 ]. Set- 
ting Eq = nuj e g then gives the roton occupation number 
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The large trap frequency enhancement rj resulting from 
the COM coupling enables one to reach the otherwise 
very elusive quantum regime n < 1. Without that effect, 
a typical temperature of T ~ 2 /iK 2l| gives n~4x 10 3 
for q/c2 = 10~ 3 , indicative of large thermal excitations 
of the rotor. The COM enhancement, with typical values 
of U ~ 2tt x 100 Hz, c 2 - 2?r x 20 Hz (for 23 Na atoms), 
N ~ 10 5 H, 7 ~ 27r x 50 kHz, and K ~ 2tt x 3 MHz, 
gives n m in < 1 for all suitable values of q. For very low 
temperatures T ~ 500 pK |22|, we have n m - m ~ 10~ 4 (see 
Fig. 1), illustrating how the COM technique opens up 
the possibility to observe a single-roton state 0] or even 
quantum jumps of a spin gas [2j-|j|. 

In summary, we have investigated the optomcchanics 
of an antiferromagnetic BEC by using a formal analogy 
with a torsional oscillator that is quadratically coupled to 
a single cavity mode. The resulting bistability-free effect, 
the noise spectra of both the atoms and the transmitted 
photons, and the quantized energy of the atoms were dis- 
cussed. That proposal provides an ideal nondestructive 
tool for the control of quantum spin dynamics , and 
facilitates the experimental study of the deep quantum 
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